Analysis of Ultrafast Transient Absorption Spectroscopy Data using Integrative Data Analysis Platform: 1. Data Processing, Fitting, and Model Selection

This manuscript describes a workflow for analysis of transient absorption (TA) spectroscopy data using Integrative Data Analysis Platforms (IDAP) software package. Time-dependent spectral series are analyzed through evaluation of the isosbestic point and kinetics of excited state and ground-state bleach decays. Model fitting and selection based on Akaike\u27s Information Criterion is discussed. As a practical example, we analyze excitation decays of a common protein label, Alexa Fluor 647

Resolving Three-State Ligand-Binding Mechanisms By Isothermal Titration Calorimetry: A Simulation Study

In this paper, I theoretically analyzed ITC profiles for three-state equilibria involving ligand binding coupled to isomerization or dimerization transitions. Simulations demonstrate that the mechanisms where the free or ligand-bound protein undergoes dimerization (such that the ligand cannot bind to or dissociate from the dimer) produce very distinctive titration profiles. In contrast, profiles of the pre-existing equilibrium or induced-fit models cannot be distinguished from a simple two-state process, requiring data from additional techniques to positively identify these mechanisms

Folding Under Inequilibrium Conditions as a Possible Reason for Partial Irreversibility of Heat-Denatured Proteins: Computer Simulation Study

Using computer simulations we have studied possible effects of heating and cooling at different scan rates on unfolding and refolding of macromolecules. We have shown that even the simplest two-state reversible transition can behave irreversibly when an unfavorable combination of cooling rate, relaxation time and activation energy of refolding occurs. On the basis of this finding we suppose that apparent irreversibility of some proteins denatured by heat may result from slow relaxation on cooling rather than thermodynamic instability and/or irreversible alterations of the polypeptide chain. Using this kinetic reversible two-state model, we estimated the effects of the scan rate and kinetic parameters of the macromolecule on its unfolding–refolding process. A few recommendations are suggested on how to reach maximal possible recovery after denaturation if refolding appears to be under kinetic control

Characterization of the Transition State of Functional Enzyme Dynamics

Through characterization of the solvent isotope effect on protein dynamics, we have examined determinants of the rate limitation to enzyme catalysis. A global conformational change in Ribonuclease A limits the overall rate of catalytic turnover. Here we show that this motion is sensitive to solvent deuterium content; the isotope effect is 2.2, a value equivalent to the isotope effect on the catalytic rate constant. We further demonstrate that the protein motion possesses a linear proton inventory plot, indicating that a single proton is transferred in the transition state. These results provide compelling evidence for close coupling between enzyme dynamics and function and demonstrate that characterization of the transition state for protein motion in atomic detail is experimentally accessible

Assignments of backbone 1H, 13C and 15N resonances in H-Ras (1–166) complexed with GppNHp at physiological pH

The small GTPase Ras is an important signaling molecule acting as a molecular switch in eukaryotic cells. Recent findings of global conformational exchange and a putative allosteric binding site in the G domain of Ras opened an avenue to understanding novel aspects of Ras function. To facilitate detailed NMR studies of Ras in physiological solution conditions, we performed backbone resonance assignments of Ras bound to slowly hydrolysable GTP mimic, guanosine 5′-[ß, γ-imido]triphosphate at pH 7.2. Out of 163 non-proline residues of the G domain, signals from backbone amide proton, nitrogen and carbon spins of 127 residues were confidently assigned with the remaining unassigned residues mostly located at the exchange-broadened effectors interface

On the Stabilizing Action of Protein Denaturants: Acetonitrile Effect on Stability of Lysozyme in Aqueous Solutions

Stability of hen lysozyme in the presence of acetonitrile (MeCN) at different pH values of the medium was studied by scanning microcalorimetry with a special emphasis on determination of reliable values of the denaturational heat capacity change. It was found that the temperature of denaturation decreases on addition of MeCN. However, the free energy extrapolation showed that below room temperature the thermodynamic stability increases at low concentrations of MeCN in spite of the general destabilizing effect at higher concentrations and temperatures. Charge-induced contribution to this stabilization was shown to be negligible (no pH-dependence was found); therefore, the most probable cause for the phenomenon is an increase of hydrophobic interactions at low temperatures in aqueous solutions containing small amounts of the organic additive. The difference in preferential solvation of native and denatured states of lysozyme was calculated from the stabilization free energy data. It was found that the change in preferential solvation strongly depends on the temperature in the water-rich region. At the higher MeCN content this dependence decreases until, at 0.06 mole fractions of MeCN, the difference in the preferential solvation between native and denatured lysozyme becomes independent of the temperature over a range of 60 K. The importance of taking into account non-ideality of a mixed solution, when analyzing preferential solvation phenomena was emphasized

Characterization of the Transition State of Functional Enzyme Dynamics

Through characterization of the solvent isotope effect on protein dynamics, we have examined determinants of the rate limitation to enzyme catalysis. A global conformational change in Ribonuclease A limits the overall rate of catalytic turnover. Here we show that this motion is sensitive to solvent deuterium content; the isotope effect is 2.2, a value equivalent to the isotope effect on the catalytic rate constant. We further demonstrate that the protein motion possesses a linear proton inventory plot, indicating that a single proton is transferred in the transition state. These results provide compelling evidence for close coupling between enzyme dynamics and function and demonstrate that characterization of the transition state for protein motion in atomic detail is experimentally accessible

NMR Studies of Escherichia Coli Acyl Carrier Protein: Dynamic and Structural Differences of the Apo- and Holo-forms

Two indicators of conformational variability of Escherichia coli acyl carrier protein (ACP) have been investigated, namely backbone dynamics and chemical shift variations of ACP. Hydrophobic interactions between the 4′-PP prosthetic group and the hydrophobic pocket enclosed by the amphipathic helices resulted in chemical shift perturbations in the residues near the prosthetic group binding sites and contact sites in the hydrophobic pockets upon conversion from apo- to holo-forms. At pH 7.9, destabilization of ACP due to negative charge repulsions and the deprotonated state of His 75 resulted in observed chemical shift changes in the C-terminal region. Model-free analysis showed that the α1α2 loop region near the prosthetic group binding site in ACP shows the greatest flexibility (lowest S2 values) and this result may suggest these flexibilities are required for structural rearrangements when the acyl chain binds to the prosthetic group of ACP. Flexibility of ACP shown in this study is essential for its ability to interact with functionally different enzyme partners specifically and weakly in the rapid delivery of acyl chain from one partner to another

The Ras G Domain Lacks the Intrinsic Propensity to Form Dimers

Ras GTPase is a molecular switch controlling a number of cellular pathways including growth, proliferation, differentiation, and apoptosis. Recent reports indicated that Ras undergoes dimerization at the membrane surface through protein-protein interactions. If firmly established this property of Ras would require profound reassessment of a large amount of published data and modification of the Ras signaling paradigm. One proposed mechanism of dimerization involves formation of salt bridges between the two GTPase domains (G domains) leading to formation of a compact dimer as observed in Ras crystal structures. In this work, we interrogated the intrinsic ability of Ras to self-associate in solution by creating conditions of high local concentration through irreversibly tethering the two G domains together at their unstructured C-terminal tails. We evaluated possible self-association in this inverted tandem conjugate via analysis of the time-domain fluorescence anisotropy and NMR chemical shift perturbations. We did not observe the increased rotational correlation time expected for the G domain dimer. Variation of the ionic strength (to modulate stability of the salt bridges) did not affect the rotational correlation time in the tandem further supporting independent rotational diffusion of two G domains. In a parallel line of experiments to detect and map weak self-association of the G domains, we analyzed NMR chemical shifts perturbations at a number of sites near the crystallographic dimer interface. The nearly complete lack of chemical shift perturbations in the tandem construct supported a simple model with the independent G domains repelled from each other by their overall negative charge. These results lead us to the conclusion that self-association of the G domains cannot be responsible for homodimerization of Ras reported in the literature

Faithful Estimation of Dynamics Parameters from CPMG Relaxation Dispersion Measurements

This work examines the robustness of fitting of parameters describing conformational exchange (kex, pa/b, and Δω) processes from CPMG relaxation dispersion data. We have analyzed the equations describing conformational exchange processes for the intrinsic inter-dependence of their parameters that leads to the existence of multiple equivalent solutions, which equally satisfy the experimental data. We have used Monte-Carlo simulations and fitting to the synthetic data sets as well as the direct 3-D mapping of the parameter space of kex, pa/b, and Δω to quantitatively assess the degree of the parameter inter-dependence. The demonstrated high correlation between parameters can preclude accurate dynamics parameter estimation from NMR spin-relaxation data obtained at a single static magnetic field. The strong parameter inter-dependence can readily be overcome through acquisition of spin-relaxation data at more than one static magnetic field thereby allowing accurate assessment of conformational exchange properties